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Abstract 

We show that under reasonably general assumptions, the first order asymptotics of the free 
energy of matrix models are generating functions for colored planar maps. This is based on the fact 
that solutions of the Schwinger-Dyson equations are, by nature, generating functions for enumer¬ 
ating planar maps, a remark which bypasses the use of Gaussian calculus. 
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1 Introduction 


It has long been used in combinatorics and physics that moments of Gaussian matrices have a 
valuable combinatorial interpretation. The first result in this direction was due to Wigner m who 
proved that the trace of even moments of a, N x N Hermitian matrix A with i.i.d centered entries 
with covariance N~^ converge as N goes to infinity towards the Gatalan numbers which enumerate 
non crossing partitions. If one restricts to Gaussian entries, that is matrices following the law /ijv 
of the GUE which is the law on the set TC]\f of A x A Hermitian matrices with density 

fiN{dA) = n diRe(Aij} dQm{Aij), 

it occurs that the corrections to this convergence count graphs which can be embedded on surface 
of higher genus, a fact which was used by Harer and Zagier m- This enumerative property was 
fully developed after’t Hooft, who considered generating functions of such moments. 

For instance, c.f Zvonkin m, we have the formal expansion 

FN{tx^) = ^log f 

k>l g>0 


with 


C{k,g) = Card{ maps with genus g 

with k stars of valence 4} 
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Here, maps are connected oriented diagrams which can be embedded into a surface of genus g 
in such a way that edges do not cross and the faces of the graph (which are defined by following 
the boundary of the graph) are homeomorphic to a disc. The counting is done up to equivalent 
classes, i.e. up to homeomorphism. Let us stress that the above equality is only formal and should 
be understood in the sense that all the derivatives at the origin on both sides of the equality match, 
it means that, for all A: G N, 


9>0 

which can be proved thanks to Wick’s formula (note above that the sum is in fact finite). 

Such expansions can be generalized to arbitrary polynomial functions (to enumerate maps with 
vertices of different degrees) and to several-matrices integrals which allow to enumerate colored 
maps; if H is a self-adjoint polynomial of m non-commutative variables, 

Tjv(H) = ^log j • • • d/XAr(Am) 

expands, when V = Y17=i some monomials Qi and real parameters ti, with q* being 

the adjoint of qi (see section El into an enumeration of colored maps. 

Our aim is to look beyond this formal work and gives a rigorous proof of this expansion. 

In the case of one matrix integrals, this problem is quite well understood at any level of the 
expansion and for any reasonable potentials V (see ^ and m for instance). 

Several matrix models are much harder. In the physics literature, the focus is mostly on a 
few specific integrals; we refer the interested reader to the reviews EH Eg. In the mathematical 
literature, fewer matrix integrals could be analyzed and only their first order asymptotics could be 
derived (see Mehta et al. El Eg and Guionnet El El). In free probability, even the problem of the 
existence of the free energy is wide open for reasonably general potential V and its solution would 
have important consequences. In combinatorics, another road was opened by Bousquet-Melou and 
Schaeffer El; following the ideas of Tutte m, by using directly bijection between maps and well 
labeled trees. 

In this paper, we shall rather focus on an even more interesting quantity than the free energy, 
namely, the limiting empirical distribution of matrices; for ^i,... , Am G it is defined as the 
linear form on the set C(Xi,... ,Xm) of polynomials of m non-commutative variables so that 

,AmiP) = ^tr(P(Ai,... ,Am)). 

Let /jy be the law on given by 

m 

i=l 


with F]sr{V) as above. 

For the following, we take a potential V = Vj = + Qi)- Then we proceed in two steps 

to relate the hrst asymptotic of gy to the enumeration of planar graphs. 
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• First, we study the solution r G C(Xi, • • • ,Xm)* of the so-called Schwinger-Dyson equations 

SD[V] : 


Tv{{Xi ViV)P) = Tv ® Tv{DiP) 

for all P G ,Xm) and i G {I,-- - ,w,}. Here, Di and Vi are respectively the non- 

commutative derivative and the cyclic derivative with respect to the variable (see paragraph 
More precisely we’re interested in properties of unicity and existence for the solution to this 
equation. 

Moreover, let us associate to branches of different colors, and to a monomial 

q(X.) = • • • Xi^ a star with p colored branches by ordering clockwise the branches corresponding 

to Xi ^, • • • , Xi^. Such a star is said to be of type q. Note that it has a distinguished branch, the first 
one, Xi^, and its branches are oriented by the above clockwise order (one should imagine the star 
to be fat, each branch made of two parallel segments which have opposite orientation, the whole 
orientation being given by the clockwise order). This defines a bijection between non-commutative 
monomial and oriented stars with colored branches and one distinguished branch, (see a precise 
description of the planar maps we enumerate in subsection 12.51) . 

We can now relate Schwinger Dyson’s equation and maps enumeration : 

Theorem 1.1 Let R > 2, then there exists an open neighborhood U C M"" of the origin (a ball of 
positive radius) such that : 

- For t G U, there exists a unique tj G C(Xi, • • • , Xm)* which is a solution to SD[l/j] and such 
that for all p, for all ii, - ■ ■ ,ip in [|1, m\], Tj{Xi^ ■ ■ ■ Wi) < 

- For all P monomial in C(JTi, • • • , Xm), t tj{P) is analytic on U and for all ki, - ■ ■ ,kn 

integers, • • • d'ff tj{P)\j^q is the number of maps with ki stars of type qi or q* and 

one of type P. 

• Then, we shall see (see section n that, under some appropriate assumptions on V, ... 
converges under py towards a solution ry to the Schwinger-Dyson equations SD[1/]. 

First if V is sufficiently large, the limit will solve a weak form of the Schwinger-Dyson equation 
(see section EU). Then we will consider convex potential V (see section EH) for which we have : 

Theorem 1.2 Let Ua be the set ofti’s for which Vj is convex, then there exists e > 0 such that for 
{ti)i<i<n £ HanH(0,e), Py_[p^] Converges to the unique solution to SD[Vj] as described in theorem 

o" 

Hence, {TVj)\t\<e are generating states for the enumeration of colored planar maps and Schwinger 
Dyson’s equations can be viewed as the generating differential equations to enumerate colored 
planar maps. This is due to the fact that the action of the derivatives Di and Vi on monomials, under 
the above bijection between stars and monomials, produces natural operations on planar maps. The 
fact that they are related with matrix integrals can be used to actually study these equations and 
obtain informations about their solutions, henceforth solving the related combinatorial problem. 
It further should give information over higher genus maps, a problem that we shall tackle in a 
forthcoming paper. 

Coming back to the free energy of matrix models, we conclude (see Theorem 13.Ill) that when the 
empirical distribution of matrices converges towards the solution to Schwinger-Dyson’s equations, 
the free energy is also a generating fnnction of the associated planar maps. In section EH we give 
another sufficient condition for this convergence to hold. More precisely, we argue that the whole 
machinery works for a general polynomial Vj with t sufficiently small if we add a well chosen cut-off. 
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Finally, we will apply these results to the study of Voiculescu’s microstates entropy and sum¬ 
marize direct applications to the enumeration of planar maps associated with some matrix models. 

The results of this paper are clearly known, at least at a subconscious level, by physicists but we 
could not find any proper reference on the subject. On a mathematical level, it is rather elementary 
and we hope it will demystify this interesting field of physics to mathematicians, or at least to 
probabilists. 


2 Schwinger Dyson’s equations and combinatorics 

2.1 Tracial states 

Let C(Xi, • • • ,Xm) be the set of polynomial functions in m self-adjoint non-commutative vari¬ 
ables. We endow C(Xi, • • • , X^) with the involution given for all z G C, alHi, • • • , ip G {1, • • • , m} 
and all p G N, by 

{zXi,---XiJ = zXi^---Xi,. 

We will say that P in C(Xi, • • • , Xm) is self-adjoint if P* = P. 

For any R > 0, completing C(Xi, • • • , X^) for the norm 


\p\\r = 


sup 


A c*-algebra 


sup ||T’(ai ; * * * ? )IU 


produces a C*-algebra C(Xi, • • • ,Xm)R = {C{Xi, ■ ■ ■ , 11.||ij, *). 

We let C(Xi, • • • , Xm)* be the set of real linear forms on C(Xi, • • • , Xm) (i.e linear forms such 
that T{a*) = r(a)), and denote C(Xi, • • • , Xm)*R the subset of C(Xi, • • • , Xm)* of continuous forms 
with respect to the norm ||.||r, i.e the topological dual of C(Xi, • • • ,Xm)R- 

We let Af”* be the set of laws of m bounded self-adjoint non-commutative variables, that is the 
subset of elements r of C(Ai, • • • , Xm)* such that 


t{PP*)>0, t{PQ)=t{QP) VP,Q G C(Xi,--- ,X^), r(/) = 1. (1) 

It is not hard to see that for any R < oo, Al^ = C(Ai, • • • , Al™ is a compact metric space 

for the weak*-topology by Banach-Alaoglu theorem. Elements of Al™' = U/{>oAl^ are said to be 
compactly supported, by analogy with the case m = 1 where they are indeed compactly supported 
probability measures. A family (rtjtg/ of elements of Af^ for some i? < oo is said to be uniformly 
compactly supported. 

To deal with variables which do not have all their moments, we eventually can change the set 
of test functions and, following [H], consider instead of C(Ai,--- ,Xm) the complex vector space 
C™(C) generated by the Stieljes functionals 

^ m 

ST^iC) = { n G C\M,«f G M,p G N} (2) 

l<i<p k=l 

where is the non-commutative product. We can give to ST*^{C) an involution and a norm 

||P||oo= sup sup ||P(ai,--- ,am)||oo 

.4C*-algebra a.i=a*&A 
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where the supremum is taken eventually on unbounded operators affiliated with A, which turns it 
into a C*-algebra. We will denote C^(M) = {G = F + F*,F e C™(C)}. We will let be the set 
of linear forms on C™(C) which satisfy (P) (but with functions of C^(C) instead of C(Xi, • • • , Xm))- 
If one equips with its weak topology, then is a compact metric space (see jHI)- 

2.2 Non-commutative derivatives 

We let Di : C{Xi, • • • , Xm) —> C(Xi, • • • , Xm)®C{Xi, • • • , Xm) be the non-commutative deriva¬ 
tive with respect to Xi given by the Leibnitz rule 

Di{PQ) = DiP X (1 (g) Q) -F (P (g) 1) X DiQ 

for any P,Qg C(Xi, • • • , Xm) and the condition 

DiXj = li=jl (g) 1 . 

In other words, if P is a non-commutative monomial 

DiP = ^ Pi (g) P 2 

P=P^XiP2 

where the sum runs over over all possible decomposition of P as PiXiP 2 . This definition can be 
extended to C^(C) by keeping the above Leibnitz rule (but with P, Q in C™(C)) and 

m mm 

Fi(^Zi ^ ^ OikX-k) — ^ ^ {Zi ^ ^ 

k=l k=l k=l 

We also define the cyclic derivative Vi as follows. Let m : C(Xi, • • • ,Xm) <8> C(Xi, • • • ,Xm) —> 
C(Xi, • • • ,Xm) (resp. C^(C) (g) C™(C) —> C^(C)) be defined by m{P ® Q) = QP. Then, we set 

Vi = m o Di. 

If P is a non-commutative monomial, we have 

ViP= Y. ^2^1. 

P=PlXiP2 


2.3 Schwinger-Dyson’s equation 

Let V be self-adjoint and consider the following equation on C(Xi,--- ,Xm)*] we say that 
r e C(Xi,--- ,Xm)* satisfies the Schwinger-Dyson equation with potential V, denoted in short 
SD[V], if and only if for all f G {1, • • • ,m} and P G C(Xi, • • • ,Xm), 

r(/) = l, T0T{DiP) = T{{ViV + X,)P) SD[V] 

These equations are called Schwinger-Dyson’s equations in physics, but in free probability, 
one would rather say that the conjugate variable (or alternatively the non-commutative Hilbert 
transform) D*1 under r is equal to Xi + ViV for all i G {1, • • • , m}. 
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2.4 Uniqueness of the solutions to Schwinger-Dyson’s equations for small pa¬ 
rameters 

In this paper, we shall restrict ourselves to non-oscillatory integrals, that is to the case where 
tr(l/ {Xi, • • • , Xm)) is real for any m-uple of Hermitian matrices. In other words, 

tr(y(Xi,--- ,Xm))=tT{V*{Xi,--- ,Xm))=tT{2-^{V + V*){X^,--- ,X^)) 

for any m-uple of Hermitian matrices. Thus, we shall assume that 

n 

V{Xi,--- = - ,X^) = Y,ti{qi{Xu--- - ,X^)) 

i=l 

where the qiS are monomial functions of m non-commutative indeterminates and t = (fi, • • • ,tn) 
are real parameters. 

In this paragraph, we shall consider solutions to SD[I^] which satisfy a compactness condition 
that we shall discuss in the following subsections. Let R G M'’“ (We will always assume R > 1 
without loss of generality). 

(H(R))4n element r G C(Xi, • • • ,Xm)* satisfies (H(R)) iff for all /c G N, 

max |t(Wi • •-XiJI < 


In the sequel, we denote D the degree of V, that is the maximal degree of the g's; qi{X) = 
X,i ■ • • X,i with, for 1 < i < n, deg(g'i) =: < D and equality holds for some i. 

1 dj 

The main result of this paragraph is 

Theorem 2.1 For all t G M"', there exists A{t) = ^(|f|) G M+ with \t\ = maxi<i<„ \ti\, A{\t\) goes 
to infinity when \t\ goes to zero, so that for R < ^(|t|), there exists at most one solution tj to 
SD[1^] which satisfies (H(R)). 

Remark : Note here that it could be believed at first sight that the solutions to SD[V] are 
not unique since they depend on the trace of high moments r(gjP). However, our compactness 
assumption (H(R)) gives uniqueness because it forces the solution to be in a small neighborhood 
of the law tq = fr™ of m free semi-circular variables, so that perturbation analysis applies. We shall 
see in Theorem lOl that this solution is actually the one which is related with the enumeration of 
maps. 

Proof. 

Let us assume we have two solutions r and rb Then, by the equation SD[V], for any monomial 
function P of degree I — 1, for i G {1, • • • , m}, 

(r - fffiXiP) = ((r - r') 0 r)(AP) + (r' ® (r - r'))(AP) - (r - fifiV^VP) 

Hence, if we let for I G N 

sup \t{P)-t'{P)\ 

monomial P of degree 1 

we get, since if P is of degree I — 1, 

1-2 

DiP = ^pl®pf_2-k 

k=0 
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where i = 1,2 are monomial of degree k or the null monomial, and 'D^V is a finite sum of 
monomials of degree smaller than D — 1, 

= max max {|T(XjP) - r'(XjP)|} 

p of degree i-i i<*<™ 


1-2 D-1 

< 2j2^k{T,T)R'^~^~^ + C\t\ 

k=0 p=0 

with a finite constant C (which depends on n only). For 7 > 0, we set 

d'riT,T') = ^7^A/(r,r'). 
l>0 

Note that under (H(R)), this sum is finite for 7 < (R) Summing the two sides of the above 
inequality times 7 ^ we arrive at 

D-l 

d^{T,T') < - -fR)-^d^{T,T') + C\t\ 

p=0 

We finally conclude that if {R, |t|) are small enough so that we can choose 7 G (0, R~^) so that 

D-l 

2^\1 - ^R)-^ + C\t\ Y < I 

p=0 

d^{T,T') = 0 and so t = t' and we have at most one solution. Taking 7 = {2R)~^ shows that this 
is possible provided 

W+c|i| 

_ 1 

so that when |t| goes to zero, we see that we need R to be smaller than A(|t|) of order |t| 0 - 2 . 

□ 


2.5 Combinatorics 

In this paragraph we describe the combinatorial objects we are considering. 

To describe the enumeration we have to deal with, let us associate a colored star to any mono¬ 
mial. We associate to each i G {1, • • • ,m} a different color. Then, we define a bijection between 
oriented branch-colored stars with a distinguished branch and non-commutative monomials as fol¬ 
lows. For any i G {1, • • • ,m}, we associate to Xi a branch of color i. We shall say that a star is 
of type q{Xi, ■ ■ ■ , Xm) = Xi^ • • • X^ if it is a star with I branches which we color clockwise; the 
first branch will be of color fi, the second of color 12 ... etc ... until the branch is colored with 
color ii. Note that this star possesses a distinguished branch, the one corresponding to X 7 , and an 
orientation, corresponding to the clockwise order. By convention, the star of type g = 1 is simply 
a point. 

A planar map is a connected graph embedded into the sphere with colored stars, each branch 
is glued with exactly one branch of the same color and the edges obtained in this way do not 
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marks the first branch 



Fig. 1 - The star of type q{X) = XjX^XfX^ 


cross each other. Hence, branches are thought as half edges. Maps are only considered up to an 
homeomorphism of the sphere. Now we will be interested in enumerating maps with a fixed set of 
stars, we define for qi a family of non necessarily distinct monomials and ki a family of integers : 

^i)) {Q 2 , ^ 2 )) • • •) {Qu, kn)) = ti{planar maps build with ki stars of type qi}. 

We denote in short Aio{P, {qi, ki), ..., (g^, kn)) = 1), (gi, ki), ..., (g„, kn))- In that set, each 

star is labeled and has a marked branch (which correspond to its first variable) so that for example 
Mq{{X'^ ,2)) = 36. Now what we are really interested in is enumerating maps with a fixed number 
of stars of type g or q* so that we define : 


n 

Wl((gi,fei),...,(gn,A:n)) = ^ Y\Cl)Mo{{qi,pi),{ql,ki - Pi), - - - ,{qn,kn),{ql,kn - Pn)) 

2=1 

l<i<n 

and M{P, (gi, ki),{qn, kn)) = M((P, 1), (qi,ki),(qn, kn))- 

This quantity enumerates the number of ways to build a map on stars of fixed types up to the 
symmetry induced on stars by the operand *. 

Due to the fact that everything is labeled, we enumerate lots of very similar objects. A way to 
avoid this problem is to look at the maps as they are enumerated by combinatoricians (see Q). 
The idea is to forget every label and to add a root which is defined as a star and a branch of this 
star. We will say that a map is rooted at a monomial of type P if its root is of type P with the 
marked branch the first one in the above construction of a star from a monomial. We can define 
for P a monomial, ki a family of integers and qi a family of (this time) distinct monomials. 

1 ^ 0 (P) (^ 1 ) ^ 1 )) ( 92 , ^ 2 ), • • •, (Qn, kn)) = t){ rooted planar maps with ki stars of type qi 

and one of type P which is the root } 


and 


V{P,{qi,ki),---,{qn,kn)) = ^ V{P,{qi,pi),{ql,ki - Pi), - - - ,{qn,kn),{qn,kn - Pn)) 

l<i<n 
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To go from these rooted maps to the previous one we only have to label each star and be careful 
about the symmetry of the stars in order to specify a branch by star. More precisely, let us define 
the degree of symmetry s{q) of a monomial q as follows. Let lo : C(Xi, • • • , X^) C{Xi, ■ ■ ■ , X^) 
be the linear function so that for alH^ G {1, • • • , m}, 1 < k < p 

ujiXi,Xi,---Xi^)=X,,---X,^X,, 

and, with = iv o define 

s{q) = tJ {0 < p < deg (q) - l\uj^{q) = q}. 

We easily see that for all monomial P, distinct monomials qi (but eventually, one of them may be 
equal to P), and integers ki : 


V{P, {qi,ki), {q 2 , k 2 ), ■■■, {qn, fc„)) 


M{P, (gi, fci), (g2, fc2), • • •, {qn, kn)) 


( 3 ) 


2.6 Graphical interpretation of Schwinger-Dyson’s equations 

We shall now make an assumption on the solutions of Schwinger-Dyson’s equation SD[l/j] when 
the parameters belong to an open convex neighborhood of the origin, namely 

(H) There exists a convex neighborhood U G M”, a finite real number R and a family {rj, t G U} 
of linear forms on C{Xi, • • • , Xm) so that for all t in U, tj is a solution of SD[l/j] which satisfies 
H{R). 

Note that up to take a smaller set U, we can assume that the conclusions of Theorem 12.11 are 

valid, i.e R < ^(|t|) for all t G ?7 since A(|t|) blows up as |t| goes to zero. 

The central result of this article is then 

Theorem 2.2 Assume that (H) is satisfied. Then 

1. For any P G C(Wi, • • • , Xm), t € U ^ origin in the sense that for all k = 

{ki,k 2 , • • • , kn) G there exists e{ki+k 2 -\ - \-kn) > 0 so that ■ ■ ■ d^firfiP) exists on = 

Ur\B{Q, e) with B{0, e) = {t G M" : \t\ < e}. We let t^{P) = ■ ■ ■ d’lfT-fiP)\i^Q. 

Then, we have for all P G C(Xi, • • • , Xm) and all i G {1, ■ ■ ■ , m}, 

n 

T^{XiP)= ^ llClWP^T^-P{DiP)+ kjT^-^i{Viq,pViq*P) (4) 

0<pj<kj j = l 
l<j<n 

where lj{i) = li=j and t^( 1) = lp=o' 

2. For any monomial P G C(Wi, • • • , Xm), any ki, - ■ ■ , G N, 

T^{P)=M{P,{qi,ki),...,{qn,kn)). 


Proof. 

• The smoothness of t comes as in the proof of Theorem 12.11 from Schwinger-Dyson’s 

equations and induction on the degree of the test polynomial function. Denote V = Vj, t = tj and 
take t = (ti,- • • ,tn),t' = {tfit2,..,t'J G U. By SD[V], 

in - + ViVfijP] = {tj - Tjf) ® Tj{DiP) + Tjf ® [tj - Tj'){DiP) + Tj>[{ViVj> - 
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By our finite moment assumption, we deduce that if P is a monomial function of degree / — 1, for 
any i G { 1 , • • • , m}, 

\tj[{X, + V,V,)P] - + V,Vj)P]\ 

1-2 

<2^ max \ri[Q] - Tj>[Q]\R^-^-’^ + ^ 

monomial of degree <fc i<i<n 

Thus we deduce that for any p G N, 

^lipXt') = max max \Tj{XiP) - Tj'{XiP)\ 

* p monomial of degree p-i 

1—2 n 

< 2 ^ Afc(Tj, ^ ^ ^ |ti|A;+cii-i(rj, Tj/) + ^ 

fc=0 i=l l<i<n 

Now, let 7 G (0,P“^) and let’s sum both sides of this inequality multiplied by 7 ^ to obtain, with 
d'firpTj^) = E«>oVA«(Tt,rj/), 

d^{n,Tj>) < 2(1 - 

n 

+ '^Hl~'^^""^d^{Tj,Tj>) + {l--iR)-^ \ti - t'i\R^-^. 

i=l l<i<n 

Since by definition Ai{tj,Tji) < 2R\ d^{Tj,Tji) is finite for 7 P < 1 we arrive at 

(1 - 27^(1 - P7)-i - Y \U\l~^""‘^)d^{Tt,Tj:) <{i- RiY^ Y \ti - ■ 

l<i<n l<i<n 

Now, for |t| small enough, we can find 7 = 7 (|t|) > 0 so that 

1 - 27^(1|ti|7“^+">0 

l< 2 <n 


and so 

Yi^Mpxt) < c{t) Y 1 ^* “ 

l>0 l<i<n 

which implies that for alH G N 


Az(rj,rj/) < C'(t)7 ' Y 

l<i<n 


so that for any monomial function P, t ^ 'R{P) is Lipschitz in Ue := U nP(0, e) for e small enough. 
Moreover, we have proved that there exists r]o{e) = 7 “^ < 00 , so that 

A Co{e)r]o{ey\t — t'\ with |t —t'| = max |tj —t(|. (5) 

^ l<i<n 

Consequently, tj is almost surely differentiable in and the derivative satisfies 

+ PPiqkP] = dt^Tj 0 TjiDiP) + tj 0 dt^Tj{DiP) ( 6 ) 
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for almost all t ^ U^. Since C(Xi,--- ,Xjn) is countable, these equalities hold simultaneously for 
all P G C(Xi, • • • , Xm) almost surely, let C/' be this subset of 11^ of full probability. 

© implies that 


max , \dt^Tj{P)\ < Co{e)'no{ey 

L<k<mp monomial of degree i 


for all t G Uy This bound in turn shows that we can redo the argument as above to see that for \t\ 
small enough, t —> dti^Tj{P) is Lipschitz. Indeed, if we set 

^i(0 = dTj) = max mp 1 ^ 4 - dt^Tj>{P)\ 

L<k<mp monomial of degree i 


we get, for t',t G C/', 


1—2 n 

Ai(0 < 2 + Co(e)|t - f\lmiey + Y, + di - 1) 

k=0 i=l 

SO that we get that by summation, for 7 < min(i?“^,r 7 o(e)“^), 

n 

(1 - 2(1 - PjrV - Y Y < 7 'Co(e)(l - - U\. 

i=l l>0 

Hence, again, we can choose ?/i(e) < 00 big enough so that there exists C'i(e) < 00 so that if e is 
small enough 

Ai{l)<C,{e)m{ey\t-f\. 

In particular, this shows that we can extend t G U'^ ^ dt^T-i{P) for all monomial functions P 
continuously in Ue and so the equality © holds everywhere. Now, we can proceed by induction 
to see that t —> tj(P) is C°° differentiable in a neighborhood of the origin. More precisely, for any 
k = {ki, - ■ ■ , kn) there exists e = e(A:i + ^2 + • • • + kn) > 0 so that on t/g, 

rf (P) = (_i)P+-+fc.afci . . . 

exists and furthermore satisfies the equation 

n _ _ 

Ty{{X,+V,Vj)P)= Y X{cyyry®Y^{D,P)+ Y k^TY^{{V,q^+V,q*)P) 

0<p^<k^ i=l i<j<rn 

l<i<n 

Applying this resnlt at the origin, we obtain the announced result. 

• We finally show the combinatorial interpretation of 0 . 

Let us first interpret graphically = tq. tq satisfies by definition SD[0] which is well known 
to have a unique solution given by the law of m free semi-circular variables (see Voiculescu EZl). 
Then, ro(Aij ■ ■ ■ Xi^) can be computed for instance using cnmulants techniques as developed by 
R. Speicher Uni ; it counts the number of planar maps which can be constructed from the star 
associated to Xi^ ■ ■ ■ Xi^ by gluing together the edges of the star of the same color. A way to prove 
that is to remark first that if we have a star with two branches of the same color, there is only 
1 = T^(l) ways to glue them. We then proceed by induction over the degree of the monomial 
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function. We let M(P) be the number of planar maps with labeled stars with a star of type P 
and shall show that it satisfies the same induction relation than t^(P). Let i G {1, • • • ,m} and 
P = XiQ. To compute Ai{XiQ), we break the edge between the distingnished branch Xj and 
the other branch of Q with which it was glued, then erasing these two branches. Since the maps 
are planar, this decomposes the planar map into two planar maps (see figure n corresponding 
respectively to the stars Qi, Q 2 for any possible choices of Qi, Q 2 so that Q = QiXiQ 2 . Hence 

M(XiQ) = MQi)MQ2)- 

Q=QiXiQ2 

Thus, if Xi{R) = tq{R) for all monomial of degree strictly smaller than P, 

M(XiQ) = ^ to{Qi)to{Q2) = To 0 To{DiQ) 

Q=QiXiQ2 

which completes the argument since the right hand side is exactly ro(XjQ). 



Edge to be broken 


Boundary of the maps 


Fig. 2 - The decomposition P{X) = XiX|XfX| into XiX^Xi (g) XfX^ 

We now consider the general case; let us assume that for |A;| < M, the graphical interpretation 
has been obtained for all monomial and that for \k\ = M 1, it has been proved for monomial 
of degree smaller are equal to L. By the preceding, we can take M > 1 and L > 1 since for all 
A: 7 ^ 0, r^(l) = 0. Again, we shall show that Xi{P, {qi, fci), • • • , {qn, kn)) satisfies the same induction 
relation than t^{P). 

Let us consider a star of type XiP (rooted at the branch Xj, with its inner orientation) 
with P a monomial of degree less than L and \k\ = Y^ki = M + 1. Now, in order to compnte 
M.[XiP, {qi, ki), • • • , {qn, kn)), we break the edge between the distinguished branch Xj (which has 
color i) and the other branch with which it was glued. 

The first possibility is that it was glued with an edge of the star P. Then, since the maps are 
planar, this decomposes the map in two planar maps. If this branch was given by the Xj so that 
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P = PiXiP 2 , one of this planar map contain the star of type Pi and the other the star of type P 2 , 
which have also a distinguished branch and are oriented. If one of this planar map is glued with 
kj stars of type Qj or q* , 0 < kj < n, the other map is glued with the remaining stars, that is 
kj — pj stars of type Qi or q*. There are 0^=1 Ck ways to choose pj among kj stars of type qj or q* 
for 1 < j < n (recall here that stars are labeled). Since we do that for all (Pi, P 2 ) so that P have 
the above decomposition, we obtain the planar maps corresponding actually to the stars associated 
with the monomials of DiP. Note that the case where one of the monomial in DiP is the monomial 
1 shows up when P = XiQ or QXi for some monomial Q and the weight corresponds then to the 
case where we glue the first branch Xi in XiP with its left or right neighbor. In this case, none of 
these two branches can be glued with another star, and there is only one possibility to glue these 
two branches otherwise, which corresponds to the weight t^( 1 ) = I^^q. 

Hence, the number of planar maps corresponding to this configuration is given by 


E E n Cl]M[Pl,{qi,Pl), - ■ ■ ,{qn.,Pn))M{P 2 ,{qi,ki -pi),--- ■,{qn,kn - Pn)) 

0 <pj<kj P=P^XiP 2 l<j<n 

l<j<n 

= E n C«rI'®r‘-i>(AP) 

l<j<n 

where we finally used our induction hypothesis. 

The other possibility is that this edge is glued with a star of type ( 7 ^ for some j e {1, • • • , n}, e e 
{.,*}. In this case, erasing the edge means that we destroy a star of type and replace the stars 
of type XiP and q glued together with a single star P glued with the star in place of Xi with an 
edge of color i removed; if g| = QiXiQ 2 , we replace the two stars of type XiP and q^- by a single 
one of type Q 2 Q 1 P (see figure | 21 ). Since we do that with all the possible edges of color i in q^ we 
find that we can glue all monomials appearing in Piqj, and so the corresponding weight is given 

by {PiqljP) times kj, the number of ways to choose one star among kj of type g^. 

Hence, by induction, we proved that the number of planar maps with kj stars of type qj or q* 
and one of type XiP is given by 


M{{XiP, 1), (gi, /ci),..., (gn, kn)) 


E E k,T~^-^^mq,+q*)P) 

0<Pj<kj j=l 

1 < j < m 


tHx,P) 


( 7 ) 

( 8 ) 


for alH G {1, • • • , m}. This shows that the graphical interpretation holds for all L and |A:| < M + 1. 
We can start the induction since we know that t^(1) = This completes the proof. 

□ 


Remarks : 

1. This graphical interpretation can be sometimes simplified for particular V. For example, 
consider V),u,c = tA^ + uB^ + cAB which appears in the Ising model. First, one may notice 
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X2 



Fig. 3 - The merging of q{X) = XfX^ = XiXiX^ and XiF into XiX^F 


that this is not in the form of the theorem but we can replace tr(F) by 

tr(^(A^ + + l(AB + BA)} 

so that the theorem can be applied. Now, one may see that there will some redundancy in the 
enumeration of maps given by this potential as for example vertices of type AB are isomorphic 
to vertices of type BA. But there will be some simplification with the factors ^ so that finally 
the theorem will give, for all monomial F, if Tt^u,c is the a solution to SD[Vt,u,c] given by the 
theorem then 

= MoiF, {A\m), (B^n), iAB,r)) 

2. Note that this graphical approach can be generalized to matrix models with more complex 
potentials involving tensor products. For example, one can consider a potential V which is a 
sum of monomials and of tensor products of monomials : 

F* = ^ tiiqj 0...0qf + {q})* 0 ... ® (qf)*) 


and the associated measure with density with respect to given by Zjy^e ‘*(tr)®‘^Vi^ 
Then one can write the generalized Swinger Dyson’s equation : 

T (g) r( A-P) = T{XiF) + {ql®---® FiqlF (g) • • • (g) 

k,3 

qt + {qiT ® ® MifP®---® {qtT) 

The previous results remain valid up to a graphical interpretation of the new term. For 
example q^ ® ® q^ will be a bunch of k loops, the first one containing the branches of q^ 
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in the clockwise order, the second one the branches of ... Note that one of these k loops is 
redundant, so that we can choose to delete one loop in the graphical representations. These 
vertices will split the map into parts as the vertices which will be placed in a loop can not be 
linked to any vertices in an other loop. 


2.7 Existence of an analytic solution to Schwinger-Dyson’s equation 

The aim of this section is to prove that for all monomials {(lj)i<j<n) there exists a convex 
neighborhood of the origin (actually an open ball) and a finite constant R so that hypothesis (H) 
of section El is satisfied. Moreover, we show it depends analytically on t in a neighborhood of the 
origin. Let Vj be as before. 

Theorem 2.3 There exists an open neighborhood U C M” of the origin (a ball of positive radius) 
such that for t G U, there exists tj G C(Xi, • • • , Xm)* satisfying SD[l/j] such that : 

- t ^ Tj is analytic on U, i.e. there exists r^, k G N"' in C(Xi, • • • , Xm)* such that for all P in 
C(Xi, • • • ,Xm), t in U, 

n 

l<2<n 

and the serie converges absolutely on U. 

- = ,{qn,kn)) 

- There exists R < oo so that for all t G U, all ii ■ ■ ■ ii G {I, ■ ■ ■ , m}*, all I G N, 

\TjiX,,---X,,)\<R^. 


Remark : It will be useful to use sometimes for an alternative expression related to the 
enumeration of rooted maps. Using ©, one can obtain inside the domain of convergence, for all 
monomial P : 

'R{P) = n ■ ■ ,{qn,kn)). 

fcgpjn l<i<n 

Proof. 

If we have such a solution, it satisfies assumption (H) and by ( 0 ) if /c! = 


T^XiP) 

k\ 


E 

0 <Pj<kj 

l<j<n 


E 

P=PlXiP2 


tP(Pi) t^-p{P2) 
P- {k — p)\ 


+ E 

1 < j < m 
kj 7^0 


r"-U(p,(g^. + g*)P) 

(k-ip\ 


where the second sum runs over all monomials Pi, P 2 so that P decomposes into PiXiP 2 . We can 
use this formula to define the by induction, the graphical interpretation is directly satished. 

We must control the growth of the t^(P)’s. Our induction hypothesis will be that for k so that 
Yli ki < M — 1 and all monomial P, as well as for '^ki = M and monomials P of degree smaller 
than L, 

< ^E kij^degP n Ck i CdegP 

i 

where the are the Catalan’s numbers which satisfy 


r^(P) 

k\ 


Ck+i ^ ^ CpCj^—p, Cq 1 , 


p=0 


V/,A:eN. 

Cl 


( 9 ) 
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Here, degP denotes the degree of the monomial P and we can assume B > 2 without loss of 
generality. Our induction is trivially true for k = 0 and all L since = cr™' is the law of m free 
semi-circular variables which are uniformly bounded by 2 so that 

|r°(P)| < 

Moreover, it is satisfied for all k and L = 0 since then r^(l) = 1 ^=q. Let us assume that it is true 
for all k such that ^ fcj < M — 1 and all monomials, and for k such that 'Y^ki = M and monomials 
P of degree less than L for some L > 0. Then 


s E E ^E kij^degP-l n CpjCk,-p, CdegPi C'degP2 

Q<Pj<kj P=PiXiP2 i=l 

l<j<n 

+2 n ^CdegP+deqqi-l 

l<l<n j 

< A^>^^B<^^sP+lYlCkXdegP+l 

i 

where we used in the last line. It is now sufficient to choose A and B such that 



l<j<n 


Qdegqj-2^degqj-2' 


A 


T^(X,P) 

k\ 


An- V fid.egqj-2Adegqj-2 

^ I _ 2 _ 

H2 + A 


< 1 


(for instance B = 2”'*'^ and A = AnB^~‘^A^~‘^) to verify the induction hypothesis works for 
polynomials of all degrees (all L’s). 

Then 

fceN" *■ 

is well defined for \t\ < {AA)~^. Moreover, for all monomial P, 


n n 

niP)\ < E IK'i'i ^)fci(4^)degP < _ 4^i.)-l(45)<ie3P. 

fcgN" i=l i=l 


SO that for small t, Tj has an uniformly bounded support. 

□ 

Hence, we see that the enumeration of planar maps could be reduced to the study of Schwinger- 
Dyson’s equations SD[V]. For instance, the asymptotics of such enumeration can be obtained by 
studying the optimal domain in which the solutions are analytic. Matrix models can be useful to 
study also the solution. For instance, we shall deduce from this approach that the solutions to 
SD[V] are tracial states (the positivity condition being unclear a priori). 
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3 Existence of tracial states solutions to Schwinger-Dyson’s equa¬ 
tions from matrix models 

We let y = be a polynomial function as before, and consider 

and fly the associated Gibbs measure 

fiy{dAiy ,dAm) = fiN{dAi) ■ ■ ■ fiN{dAm). 

This is well defined provided that we assume that the monomials of highest degree in Vj are 
sufficiently large to make Zy finite. We shall assume for instance that 

Fi(x)= ^ t,(g,(X) + g:(X))+ ^ uxp_^ (10) 

l<i<n n-\-l<i<n-\-m 

with D even and monomial functions qi of degree less or equal than D — 1 and ti > 0 for i £ 
{n + 1, • • • , n + m}. We shall see in the last paragraph of this section that such assumption can be 
removed provided a cut-off is added. 

The empirical distribution of m matrices A = (Ai, • • • , A^) £ is dehned as the element of 
such that 

=/iM(i") = ^tr(T(Ai,... ,Am)) 

for all F £ C^(C). Note that the empirical distribution could be defined as well as an element of 
Al"* but since the random matrices (Ai, • • • , Am) under fiy have a priori no uniformly bounded 
spectral radius, the topology of weak convergence would not be suitable then. 

We shall see that if we know that a limit point of under fiy are compactly supported, then it 
satishes SD[V]. In a second part, we shall give examples of potential V for which this assumption is 
satished. Finally, we discuss localized matrix integrals and show that bounded solutions to SD[V] 
for small potentials can always be constructed by localized matrix integrals. 

3.1 Limit points of empirical distribution of matrices following matrix models 
satisfy the SD[V] equations 

We claim that 

Theorem 3.1 Assume (jini) . Then 

1. There exists M < oo so that 

limsup < M 

N^OO 

fly almost surely for all i £ {1, ■ ■ ■ , m}. 

2. The limit points of ... for the {C) -topology satisfy the ‘weak’ Schwinger-Dyson equa¬ 
tion 

T ® T{DiF) = T{{ViV + Xi)F) (WSD)[V] 

for all F £ Cft{C). 
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Note here that {ViV + Xi)F does not belong to C^(C) so that it is not clear what (WSD)[V] 
means a priori. We define it by the following; there exists a sequence £ C^(C) so that 

lim max sup riWiV^ — ViV — XA) = Q 
5^0 l<i<m r(xp)<M 

from which, since any F £ C™(C) is uniformly bounded, 

lim max sup \T{FViV^) — T{F{ViV + Xj))| = 0 


is well defined. 

Proof. 

• The first point is trivial since by Jensen’s inequality, 

Zn > exp{-A^^ f ^tr(F(A)) d^iNiA)} > exp{cA^} 

l<i<m 

for some c > — oo, where the last inequality comes from the fact that (see 123) 

n dnN{Ai) = a"^{V) < oo 

l<i<m 

where is the law of m free semi-circular variables. 

Now, observe that by Holder’s inequality, 

\f^AiQi)\ < m.ax fi’l{\Xi\^-^ + 1) 


so that we deduce 

m 

i^i+nf^Ai^F) - c{t)fiAi\Ai\^~F - C{t)) 

i=l 

with a finite constant c{t). Since > 0, we conclude that fi^{V) > m\t\M/2 when 
maxi<j<m ^ ^ large enough. Thus 

I.?, (,mjx f.l(Af’) > m) < (11) 

\^l<i<m J 

goes to zero exponentially fast when M > The claim follows by Borel-Cantelli’s lemma. 

• We proceed as in [^, following a common idea in physics, which is to make, in Zy , the change 
of variables Xi —>■ Xi + N~^F(K) for a given i £ {1, • • • ,m} and F £ C^(M). Noticing that the 
Jacobian for this change of variable is 

|j| ^ 


we get that 
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from which we deduce that 


/ g(7V/i^®A^(r>iF)-7V2A^(Af-iXiF{X)+l/(Xi+7V-iF(X))-y(XO)^7V^^^^ ^ 

J maxi<i<^ Aa (^f )<^ 

Hence, we conclude by Chebychev inequality and m that for M big enough, any 5 > 0, there 
exists 7? > 0, so that if we denote 

^iv = Aa ® Aa(AF) - Nfl^{N-^XiF{X) + V{Xi + iV-iF(X)) - V{Xi)) 

then 

<M}n {I^ItvI < <5}^ > 1 - e-^^. 

Moreover, 

+ N-^F[X)) - V{X,))) = + Rn 

with a rest of order X~‘^ maxi<j<m fJ-^{X[’~‘^) which we can neglect on 
maxi<j<m Aa(^?) — This shows, by Borel-Cantelli’s Lemma, that for all F G C^(M), 

Aa ® Ai(AT) - Ai(AT + ViVF) 

goes to zero almost surely. This result extends to T G C'^{C) since it can always be decomposed 
into the sum of two elements of C^(M). Moreover, if we let Aj = Hj(l + eAf)~^ = Ai(y/—1 + 
y/eAi)~^{—^/—l+y/eAi)~^ G C^(C), then again by Holder’s inequality r(|PiH(Mj)—AT(H|)|) goes 
to zero uniformly on maxi<j<m T{Af) < M. This shows that fj, —> ^{{T>iV+Xi)F) is continuous for 
the weak C^(C)-topology on {^{Af ) < M} for any F G C^(C). Therefore, since is compact, 

we conclude that any limit point of Aa satisfies 

T®T{DiF) = T{{X, + ViV)F) 


□ 


We therefore have the 

Corollary 3.2 Assume that there exists a limit point ry of Aa under Hy which is compactly 
supported. Then, it satisfies Schwinger-Dyson’s equation SD[V]. 

Proof. 

The proof is straightforward since if ry is compactly supported it is equivalent to say that Ty 
satisfies WSD[V] or SD[V] since C^(C) is dense in the set of polynomial functions (approximate 
the Afs by the Af^s defined in the previous proof). 

□ 

Let us also give the final argument to deduce convergence of the free energy from the previous 
considerations. 


Theorem 3.3 1. Assume that fi^ converges in almost surely or in expectation under 

fj,y_ towards tj solution to SD[kj] for t in a convex neighborhood U of the origin. Assume 


furthermore that maxp jXy 




(Aa(|Ap|^)) is uniformly bounded for 
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(ti,--- s,tk+i, ■ ■ ■ ,tn) £ U and N large enough for some I strictly greater than the 

degree of Vj. Then, 


F^,, = iV-2log 


7N 


'V, 


7N 

^Vj 




converges as N goes to infinity towards a limit ■ Moreover, 


r^k 

Fvj,k = - + ,tr,)i9k + qk)ds. 

Jo 


If furthermore tj is uniformly compactly supported in U, we deduce that t —> Fy- is C“ in a 
neighborhood of the origin and (—• • • d^"Fy-^ k\t=o 'i'S the number of planar maps with 
Pi stars of type qi or q* when Pk > I- 

2. Assume that for t in a open convex neighborhood U of the origin the limit points of fi^ 
under ^y_ are uniformly compactly supported. Assume further that ma.Xp\gLy_{jl^{\Xpf))\ is 
uniformly bounded (independently oft^U) for N large enough and some I strictly larger 
than the degree ofVj. Then, fi^ converges p,y_-almost surely towards tj described in Theorem 
\2.A for t £ U n 5(0, e) for some e > 0 small enough and for t £ [/ n 5(0, e), 

F{^_=N-Hog{Z^_) 

converges as N goes to infinity towards 

Py-t= Y. n AimK)). 

fceN"\(o,..,o) *’ 

Note above that the last serie as a positive radius of convergence according to Theorems 12.21 
and [Q This emphasizes that the possible divergence of Fy_ does not survive the large N limit. 

Proof. 

• By differentiating N~‘^ log Zy_ with respect to tk we obtain that 

dif^N ‘^\ogZy_ =—p,y_{fL^{qk + qk))- 

But, under assumption, {fL^{qk + 9fc))AfeN converges almost surely and is uniformly integrable so 
that I^y.{fi^{qk+qk)) ^ uniformly bounded sequence which converges as N goes to infinity towards 
Tj{qk + q^.) for t € U. Integrating with respect to tk yields the convergence with Fy- as above by 
dominated convergence theorem. The last part of the first point theorem is a direct consequence of 
Theorem 12.21 

• By Corollarv fI.2l and Theorem l2.ll we see that our hypothesis implies that for t £ 5n5(0, e) for 
some e > 0, the limit points of //^ are unique and given by tj. Hence, fi^ converges in Af almost 
surely towards tj. Since we assumed our family uniformly integrable, we deduce that IJ-y-{fi^{qi+qi)) 
converges as N goes to infinity towards Tj{qi + q() for alH £ {1, • • • ,n} and we see as above that 
for alH £ {1, • • • , n — 1}, 

N 2 1 z(j I 
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converges as N goes to infinity towards a limit 


r^i 

Fvj,i = - + ^i)ds. 

Jo 


Hence, since we know that = 1, 


1 ^ \ 
lim ^logZ^ = — 

i=l 

n 

_ / 




N—>cx) N'^ 


7N 






rti 

~'^Jq no,- as,u+i,- + gi)ds 
*=1 fci,..,A:„eN"-“ *+l<i<". ^ * 


where we used in the last line Theorem 12.hi Noting that no,..,o,ki,-" ,kn) _|_ g*^ is the number of 
planar maps with kj stars of type qj or q* for j > i +1 and fej +1 stars of type qi or q* , we conclude 
the proof. 

□ 

We shall in the next section provide a generic example where the assumption of the second 
point of Theorem ESI is satisfied (in fact, a slightly different version since we do not prove that the 
almost sure limit points of satisfy our compactness assumption, but their average do, which 
still guarantees the result). 

3.2 Convex interaction models 

Let us assume that we consider a matrix model with potential V such that 

m 

C ; (Akim e (K"’)” ,A„,)) - (12) 

^ k=l 

is convex in all dimensions for some a < 1, i.e the Hessian of is non negative for all N £ N. 
An example is V of the form 

n m 

V{Ay- ,Arn)=Y,n{Y.oJ^AkfP^ 

i=l k=l 

with non-negative tj’s, integers pi’s and real a’s. Indeed, by Klein’s lemma (c.f. Cl) , since x 
is convex, 

A ^ tii^aiAifP* 

is also convex (Here A, by an abuse of notations, denotes the entries of the m-uple of matrices 
A = (Ai, • • • , A^ffi)). 

Then, we shall prove that 

Theorem 3.4 Let V be a self-adjoint polynomial function which satisfies \liAl . Then 
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There exists Ry < oo so that 


N^OO 

for all n G N and i G {1, •'' j^}- Here, Ry is uniformly bounded by some Rm when the 
quantities (a, y( 0 , 0 , 0 ), (Pjy( 0 , 0 , 0 ))i<j<m) are bounded by M. 

- /iy [Aa] *'5 tight and its limit points satisfy SD[V]. 

- Take V = Vj = be the set of U’s for which Vj satisfies UTA) for a given 

a < 1. Fore > 0 small enough, when (ti)i<j<n G UaFB{{),e), converges to the unique 

solution to SD[V]. 

- Assume that Ua contains Ui<j<n{(0, ..,0, ij, ..,tn),0 < F < 5} for 6 small enough. Then, for 
e > 0 small enough, for t G U D B{0, e), 

F{[ = N-^\og{Z(l) 

converges as N goes to infinity towards 

Hyj= Y1 n ^ M{{qi,ki)y ,{qn,kn)). 
feN"\(0,..,0) *' 

Remark : Observe that our hypothesis is verified for all quadratic interaction models such as the 
Ising model, the ( 7 -Potts model ... etc ... as soon as the self potential of each matrix is convex. 

Proof. 

We can assume without loss of generality that a = 0 since otherwise we just make a shift on 
the covariance of the matrices under The idea is to use Brascamp-Lieb inequality (c.f Harge 
m for recent improvements) which shows that since 

/(A) = Am) 

is log-concave, for all convex function g on 

= J J (13) 

with 



Here A denotes the set of entries of the matrices (Ai, • • • , A^)- Let us apply (d with fi'(A) = 
tr(A^^) which is convex by Klein’s lemma. Hence, 

(tr((Afc -E[Afc])2P)) < M7v(tr(Af)) 

where E,[Ak]{ij) = iJ,y{Ak{ij)) for 1 < i,j < N Since the right hand side is bounded by 4^ as A 
goes to infinity, we conclude that 

limsup^y [^tr((Afc - < 4?. (14) 

N^oo A 
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We now control lE[Afc] uniformly. Since the law of ^4^ is invariant by the action of the unitary group, 
we deduce that for all unitary matrix U, 

E[^fc] = E[UAkU*] = UE[Ak]U* => E[Afc] = (15) 

We now bound ^y{fi^{Xk)) independently of N. Since V is convex, there are real numbers 
(7*)i<i<m and c > —oo, 7 i = 'DiV{0, ■ ■ ■ ,0) and c = 1^(0, • • • , 0) so that for all G N and all 
matrices (^i, • • • , Am) G 


tr(y(Ai, • • • ,Am)) > tr(^ 7 jAj + c). 

i=l 

By Jensen’s inequality, we know that Z'^ > for some d < oo and so Chebychev’s inequality 

implies that for all y > 0, all A > 0, 


d'v (lAA(^fc)l > y) 


< 


< 


i=i 


+ f e-A^Er=i 7dr(A7-iVAtr(A,) ^ dfiN{Ai)] 

i=i 


Optimizing with respect to A shows that there exists ^ < oo so that 

and so 

^^vi\y'AiXk)\) = [ yv {\yAiXk)\ >y)dy< 4VA+ f e~^^y^~'^^'>dy < sVA 
J Jy>A-/A 

where we assumed N large enough in the last line. Hence, we have proved that 

limsup |/iy(/iA(^fc))| < 8\/H. (16) 

N 

Plugging this result in m and we obtain for all p > 1 : 

limsup/iy[/iA((-4fc)^^)] < 2‘^P~^limsupny[^tT{{Ak - yvi^k])'^^)] 

N^oc N^oo J'' 

+2^^’“Mimsup(py(^tr[Hfc])2P) 

N^oo JV 

< 2‘^p-Hp + 22P-1(8\/I)2P < R^P 

with Rv = 4(1 + Sx/A). To prove the convergence of remember that is tight for 

the C^(C)-topology. To study its limit point, recall f xe~^ ^'^f{x)dx = f f'(x)e~^ I’^dx so that, for 
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p E c,7(c) 


I ^tT{AkP)d^^l{A) 


[ ^A,iij)iPe 

ij 

2j^ ^ / ( ^ ‘^QiiRjj 

ij \P=QXkR 
n N \ 

-"E E ‘.E 2PjiQhjRih I d^y{A.) 

l=l qi=QXuR h=l J 

j (^(tr ® tr)(E)fcP) - ^tr(pfcyp)^ (A) 


which yields 

I + VkV)P) -fil® filiDkP)) d/x]^(A) = 0 

Now, by convexity of V we have concentration of under (since log-Sobolev inequality is sat¬ 
isfied uniformly, according to Bakry-Emery criterion, and that Herbst’s argument therefore applies, 
see |2], sections 6 and 7) : for all Lipschitz function / on the entries 




where ||/||£ is the Lipschitz constant of /. Since for P E C™(C), A ^ id^{P) is Lipschitz with 
constant of order N~^ (see ^), we conclude that since DiP E C™(C) (g) C^(C), for all P E C™(C), 


lim 
A^—>00 


Aa ® idl{DkP)d^XN{A) - /i]C[AA] ® 7t]v[AA]Pfc7’) 


= 0 . 


Thus 

limsup (^^(//^((Afc + VkV)P)) - ®/^{/[AaKA-P)) = 0 

W—>oo 

If T is a limit point of //^[Aa] weak C™(C)-topology, we can use the previous moment 

estimates to show that even though + VkV is a polynomial function, ^y{fi^{{Xk + VkV'jP)) 
converges along subsequences towards T{{Xk + T)kV)P)), and of course /^{/[AaI ® My [AAK-^fc-^) 
converges towards T®T{DkP)- Hence, we get that the limit points of /^y [AaI satisfy the WSD[V]. 
By the previous moment estimate, this limit points are compactly supported, hence they satisfy 
SD[V]. 

When V = Vj, observe that Rj is uniformly bounded when |t| < M since Vj(0, • • • ,0) and 
('DiVj;(0, • • • , 0))i<j<m depends continuously on t. Thus, the first point of the Theorem shows that 
the limit points of Mi^IAaI uniformly compactly supported. Hence, since also we have seen that 
they satisfy SD[kj], for t small enough, converges in expectation (and therefore almost surely 
by concentration), to the unique solution to SD[I^]. The last point is now a direct consequence of 
Theorem ESI 

□ 
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Hence, we see here that convex potentials have uniformly compactly supported limit distribu¬ 
tions so that we can apply the whole machinery. We strongly believe that this property extends to 
much more general potentials. However, we shall see in the next section that we can localize the 
integral to make sure that all limit points are uniformly compactly supported and still keep the 
enumerative property, hence bypassing the issue of compactness. 

3.3 The uses of diverging integrals 

In the domain of matrix models, diverging integrals are often considered. For instance, if one 
wants to consider triangulations, one would like to study the integral 

ZNitx^) = j 

which is clearly infinite if t is real. The same kind of problem arises in many other models (c.f. 
the dually weighted graph model EU). However, we shall see below that at least as far as planar 
maps are concerned, we can localize the integrals to make sense of it, while keeping its enumerative 
property. Namely, let Vj = ^ tiQi and let us consider the localized matrix integrals given, for L < oo, 
by 

‘ A\M\oo<L 

and the associated Gibbs measure 

Here, ||A||oo = maxi<j<m ||Aj||oo and ||Aj||oo denotes the spectral radius of the matrix A,. 

We shall prove that 

Theorem 3.5 There exists eo > 0 so that for e < eo, there exists Lo{e) and T(e), L(e) going to 
infinity and To(e) going to 2 as e goes to zero, so that fort G B{0,e), for all L G [Lo(e), L(e)], 

log ^ ,{qn,kn)) 

fceN"\(0,..,0) 

Moreover, under eonverges almost surely towards tj described in Theorem A2.,‘A 

This shows that, up to localization, the first order asymptotics of matrix models give the right 
enumeration for any polynomials. The diverging integrals often considered in physics should be 
therefore thought to be conveniently localized to keep their combinatorial virtue, and are then as 
good as others. In view of LemmaEini this localization procedure should not damage the rest of the 
large N expansion neither. Note that when m = 1, the localization amounts to restrict the integral 
in the domain of strict convexity of the potential, henceforth again avoiding all issues of escaping 
eigenvalues. 

Proof. 

The proof is very close to that of Theorem 13.11 except that we have to be careful to make 
perturbations which do not change the constraint ||A||oo < L. Let i G {I,-- - ,m} and consider 
the perturbation Aj — > Aj -|- N~^h{Ai) and Aj — > Aj for j i with a compactly supported 
function h which vanishes on [—R,Rfi. Then for L > R, for sufficiently large N, and ||Aj||oo < L, 
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\\Ai + N ^h{Ai)\\^ < L so that we see that the limit points of under the localized Gibbs 
measure 

• • • , Ara) = • • • dHN{Am) 

satisfy for i G {1, • • • , m}, 


/i ® fi{DMXi)) = KiViV + Xi)h{Xi)). (17) 

These limit points are also laws of operators bounded by L, but we shall see that in fact this bound 
can be improved to become independent of L for good L’s. We fix a limit point // below ; ^ \s a. 
tracial state. We proceed as before by taking h{x) = P{(pji{x)) with a polynomial test function P 
and (j)fi a smooth cutoff function with uniformly bounded derivative (say by 2). Then, 

H (g) ^{Dih{Xi)) = fi{DiP{(j)ji{Xi)) X (p'niXi) (g) 1) 

so that if P{x) = 


2k 


p=0 


Hence, we get, for i G {1, • • • , m}, 

2k n 

^(w<(>r(a)'"+') < 2^^(i<()^(w)n/u(i<(>«(A)p''-n + iti j;i^(iAg,(x)ii<(>ii(A)|2"+') 

p=0 j=0 

Note that we can bound above the last term by Holder’s inequality so that 

n 

^f,(\nq,(X)\\cPn(X.)\^'^+^) < C max {^(|X,f^+^-1),^(|X,f^)} 
j=0 


where we assumed \(t)R{x)\ < |x|. Taking (j)R{x) = x when |x| < R/2 and 4)r{x) = 0 when |x| > R, 
(pR linear in between, we can now use monotone convergence theorem (letting R going to infinity) 
to obtain 


2k 


max < 2^ max ^(|Xj|P) max p) + C\t\ max {n{X^^) + ^)) 

I < -1 < m ^ 1*^ 7 "^TTl I < T < m I < -i < m •' 




p=0 


l<j<m 




Noting that ;u(|Xj|^'''^^ ^) < ^fj,{\Xj\‘^^) since under ^ the operators are uniformly bounded 

by L, we can improve this uniform bound as follows. We make the induction hypothesis that 

Bk := max fi{\Xj\^) < CkA^ 


for k < 2r, with Ck the Catalan numbers and some R> 1. Then, by Holder’s inequality, 

2r + 2 

{B2r+l)~^ < B2r+2 < 2C'2.+iH 2’- + C|t|(l + L^-^)A^^C2r. 


We can assume without loss of generality that B 2 r+i > 1 so that we get 

B2r+2 V B2r+1 < C 2 r+lA^''^^ {2A-^ + C\t\L^) 
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and so we can continue our induction if L is not too large (L at most of order (2C|t|)“n) so that 
there is ^4 > 0 so that 2A~^ + C\t\L^ < 1. Hence all limit points are non-commutative laws of 
operators with norms uniformly bounded by 4H <C L. Thus if i? > 4a we conclude by mi) that all 
limit points satisfy SD[4^]. Moreover, Theorem 12.11 and 12.dl apply to show that there is a unique 
solution to this equation which are laws of operators bounded by R (provided R, and so L, is neither 
too big (L < 0{\t\~^ ), for uniqueness), nor too small (L > Lo{\t\) > Lq{0) = 2 for existence. 

Note here that Lo(|t|) is the smallest real number such that \tj{P)\ < for all monomial 

P, which converges to 2 as t goes to zero by Theorem 12.41 Hence, has a unique limit point, Tj, 
and thus converges towards it. 

The formula of the free energy is then derived as in Theorem ESI since L is fixed independently 
of t small enough. 

□ 

Let us remark that if we define, following Voiculescu m. a microstates r(^, n, N, e) for /i e 
n G N, iV G N, e > 0, as the set of matrices ^i,.., of Tf™ such that 


- trAr(Aji..Ajp)| < e 

for any 1 < p < n, R, .., ip G {1,.., m}^, then we have 
Lemma 3.6 For L big enough, V = Vj and t small enough, 


(18) 


lim —^ log 

TV^oo A2 ^ 


|A||„o<L 


e dpN{Ai) ■ ■ ■ dpN{Am) 


= lim lim —^ log 

e—>0,n^oo Af—>oo N‘^ 


1 


T(Tv,n,Ar,e)n||A||oo<i 


e ^^'^^^^^^^dpN{Ai)---dpN{Ar. 


= lim lim —^ log 


e ^^^^^^^^^dpN{Ai)---dpN{Am) 


e^0,n^oc N^oc N'^ Jv{Tv,n,N,e) 

Proof. 

The first equality is a direct consequence of the previous Theorem since it is equivalent to the 
fact that py2^{T{Tv,n,N,e)) goes to one. The second comes from the fact that for n greater than 
the degree of V, 


lim lim —- log 

e—»'0,n—^■oo N^oo N‘^ 


r(rv,n,A'',e)n||A||oo<i 

1 




= -Tv(y)+ lim lim ^log^^”" (r(ry,n, A, e) n ||A||oo < L) 


e—>0,n—>cxD W—>oo 


1 


= -Tv{V)+ lim lim ^log/r^”'(r(ry,n,A,e)) 

e^0,n^oci N^oc ^ 

where we used in the last equality the result of |S], which hold when ry is the law of bounded 
operators with norms strictly smaller than L. 

□ 
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Therefore, the localization should not affect the full expansion of the integral since second order 
asymptotics are usually obtained first by a localization on a microstates in order to use precise 
Laplace method’s. 

As a corollary, we also deduce that for all Vj with t small enough, the limits of empirical 
distributions of matrices given by localized matrix models provide solutions of SD[1^]. Since these 
limits have to be tracial states, we deduce that 

Corollary 3.7 The solutions compactly supported o/SD[Vj] are tracial states whent is sufficiently 
small. 

Note that if {Pi)i<i<m is the conjugate variable of a tracial state, Voiculescu m have shown that 
Pi = T>iP for 1 < i < m and some polynomial P. This fact should be compared with our graphical 
interpretation which works only because Pi is a cyclic derivative. 

4 Applications to free entropy 

Let us recall that Voiculescu’s microstates entropy is defined, for r G UnAiff, by 
X(r) = lim limsup^log/i^™ (r(r,n,iV,e) n ||A||oo < L) 

L—*oo 

with r(T, n, e, A) the microstates defined in (I18j) . Note that the original definition of Voiculescu is 
not with respect to the Gaussian measure, but with respect to the Lebesgue measure. However, 
both definitions only differ by a quadratic term (see |H]). It is an (important) open problem whether 
in general one can replace the limsup by a liminf in the definition of y. However, from the previous 
considerations, we can see the following 

Theorem 4.1 Let n G N and {qi)i<i<n be monomials in m non-commutative variables X = 
(Xi,--- ,Xm). Let Vj(X) = Theorem \2.tA we know that there exists 

e > 0 so that for |t| < e, there exists a unique solution tj to SD[Vj]. Then, also for |t| < e, 

xiTt) = lim liminf ^log^^”" (r(T, re, V,e) n {||A||oo < L}). 

e^O.n^oo N—tOO jV ^ 

L —»-oo 

Moreover, 

^ ( + ^ \ 

X(rt) = - n fc'i X] ~ M ki), • • • , (qn, kn)). 

fceN"\(o,-,o)i=i *■ \i=i / 

remark : In particular, we see as could be expected that xiTt) > solutions 

to Schwinger-Dyson’s equations for small parameters are laws of von Neumann algebras which are 
isomorphic to the free group factor. 

Proof. 

In fact, 

X{rj) = lim limsup At log [ giVtr(y,(A))-iVtr(y,(A))^^®m(A) 

N Jr(rj,n,iV,.)n{||A||oo<L} 

= tt{Vj)+ lim lim sup At log /" 

N Jr(T^,n,N,.)n{\\A\\^m 

< rj{V,) + Fi 
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with 


Fj= lim limsup^log / 

L—>oo JV—>c» y||A||oo<L 


On the other hand, 


lim lim inf ^ log /rl”" (h (r*, n, iV, e) 0 

-0,n^oo Af—>00 iV^ 


<L) 


= lim inf log /" 

t' E^O.n^oo AT^OO J-p 


r(Tj,n,A,e)n{||A||„o<i^ 


g-Atr(y,(A))^^®n^(^) 


= 'rtiVi) + Fj+ lim lim inf — log(r(rj, n, iV, e)) 
e^O.n^oo Af^OO iV ^ * 


where we used in the last term Theorem 13.51 which implies 

lim (r(rj, n, iV, e)) = 1 

N—^OO * 

for all e > 0,n E N and L large enough. Thus, we see that x is equal to its liminf definition and 
moreover 

xin) = TiiVj) + Fj. 

Now, by Theorems 13.51 and ESI 

fceN"\(o,..,o) *’ 

whereas 

= X] n /, •M((gyfci)y ,(gi-i,fci-i),(gi,fci + l),(gi+i,fcj+i),--- ,iqn,kn)) 

i=l kjSN, l<j<n 

from which the formula for xi'^t) is easily derived. 

□ 


5 Applications to the combinatorics of planar maps 

For the sake of completeness, we summarize in this last section, the results of a few papers 
devoted to the enumeration of planar maps, either by a combinatorial approach or by a matrix 
model approach. 

5.1 The one matrix case 

We now consider the simpler case m = 1 where we only have one matrix. Let VjiA) = Y^=i 
with t 2 D > 0 a polynomial potential with an even leading power. Then it has been proven in 
Theorem 5.2 that the empirical measure satisfies a large deviation principle : 
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Theorem 5.1 Let 


J{l) = J (^Y + ~ J J^og\x- y\dfj,{y)dn{x) 


and 

I{n) = J{y)- inf J(i^) 

i^eP(R) 

then the sequence of empirical measure fi^ satisfies a large deviation principle in the scale with 
good rate function I. Moreover, the minimum of I is reached at a unique probability measure pLj so 
that 

x^ f 

Y + / log I 2 / “ x\dyj{y) = Cj, p.ja.s. 

with a finite constant Cj, and where the left hand side dominates the right hand side on the whole 
real line. 

One can notice that differentiating in x the last equation, we recover the Schwinger Dyson’s 
equation. It is not sufficient in general to determine the solution uniquely; one need the inequality 
on the whole real line to fix the support of the solution. 

These questions have also been investigated with the method of orthogonal polynomials which 
give a rather sharp description of the limit measure and emphasizes a structure similar to the 
semi-circular law. More precisely Theorem 3.1 in m gives : 

Theorem 5.2 There exists t > 0 and 7 > 0 such that if for all i, \ti\ < t and t 2 D > 
then is absolutely continuous with density of the form : 

= ^l[a,b](a;)\/(a; - a){x - b)h{x) 


with 


h{z) = f 
JCi 


V^{s) 


ds 


lc{z,R) S - ^ 

where R is such that a,b G C{z,R). Besides, the boundaries a and b can be find by the equations : 



a/(s - a){b - s) 


ds = 0 



sV^js) 


ds = 2 tt 


We now look at combinatorics of the Schwinger-Dyson’s equation with one variable, for Vfix) = 
Y^=i tix'". Remember that from Theorem 12.31 pLj can be seen as the generating function of graphs 
counted by the numbers of stars of valence i : 


ZJJ / x/j. 

>fc2DSN i=l * 

Hence, Theorem 15.21 allows to estimate the numbers of one colour planar maps. A more direct 
combinatorial approach can be developed by considering for instance the dual of those graphs. The 
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dual of a graph is simply obtained by replacing each face by a star and each edge by a transverse 
edge which link the two stars which come from the face adjacent to the edge. In that operation 
each star is replaced by a face of the same valence. As we work on the sphere we can decide that 
the face which comes from the star is the external face. 

So is also the generating function of connected planar graphs with an external face of 

valence p + 1 and enumerated by the number of faces of a given valence. Those objects are classical 
ones in combinatorics and we can follow m to find an equation on these generating functions. 
The idea is to try to cut the first edge of the external face, then two cases may occur : either the 
graph is disconnected and we obtain two graphs or it isn’t disconnected and the external face has 
grown. This two cases corresponds in the dual graph to the fact that the first branch of the root is 
a loop or not which is exactly what we use to build our combinatorial interpretation so that we can 
retrieve the Schwinger Dyson’s equation from this fact. Just by using the equation given by this 
decomposition and some algebraic tools combinatoricians have solved some models. For example 
[H] gives an equation on the generating function M {u, v) of maps whose internal faces have degree 
living in a fixed set D C N and enumerated by their number of edges and the degree of the external 
face. To translate this in our framework, one can consider for a finite T) with an even maximal 
element, 

V-,{X) = Y,tdX<^ 

d&V 

Then under this potential, for small t, the limit measure pj will satisfy our combinatorial interpre¬ 
tation. Then 

M{u,v)=^fi d {XP)vP 

Now Theorem 1 of [3] states : 

Theorem 5.3 For a serie F{z) = we will note [z^]F{z) the eoeffieient a*. Then there 

exists a unique power serie R satisfying 

R = 1 — iRiv — AR2V^ 

with 

Ri = - and R 2 = - Y,[v^]{R'^) + u- 3RI 

^ idD ^ ieD 

The number nin of maps with n edges such that every degree of internal face lies in D is then 

{R2{u) + Ri{u)^){R2{u) + 9i?i(u)^) 

(n -|- l)u‘^ 

The techniques to prove these results are most often purely algebraic. The main difference in 
nature than we could meet between the approaches by matrix models or by combinatorics to these 
enumerations is that the first provides for free additional structure; it shows that these enumerations 
can be expressed in terms of a probability measure. This point generalizes to any number of colors 
where the enumeration can be expressed in terms of tracial states. One may hope that this positivity 
condition could help in solving this combinatorics problems. 


rUn = [u 
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5.2 Ising model on random graphs 

This model is defined by m = 2 and 

V{A,B) = VIs^ngiA,B) = -cAB + Vi{A) + V 2 {B). 

In the sequel, we denote in short A for Xi and B for X 2 . It is clear that for |c| < I, 1/ is a 
convex potential as defined in m if Vi,V 2 are convex (write —2AB = {A — B)'^ — A^ — B^ or 
2AB = (A + BY — A^ — B^ to see that up to a quadratic term + 2“^|c|i?^, V is convex) 

Hence we deduce from Theorem m that 

Corollary 5.4 For c G M and Vi{x) = XjLi > * = 1) 2, set Vj^^{A, B) = —cAB+Vi{A)+V 2 {B). 
Let, for (5 > 0, 1/5 = HijlO < t*- < 5} H {|c| < 1 — <5}. Then, for <5 > 0 small enough and (t, c) G Us, 
/^v(Aa) converges towards the solution o/SD[I^ J as N goes to infinity. Moreover 

kiifflD ij j' 

r^N 

and 

F{t,c)-F{t,0)= ^ Ylt^^Mo{{A^^,k])i<j<D,{B^^,k])i<j<D,{AB,r)). 

i,j j' 

r>l 

Remark : Note that we took potentials Vi and V 2 as polynomials with even powers to guarantee 
our convexity relation but this condition could easily been relaxed by taking more sophisticated 
domains than Us in which the polynomials would remain convex. 

Proof. 

This result is a consequence of Theorem I2.2L 15.41 and 15.51 Note here that the control on 
//y (iV“^tr(Ai?)) assumed in Theorem 15.51 is satisfied due to Theorem 15.41 which provides a uni¬ 
form bound when |c| < ^ for ^ < 1. 

□ 

According to the graphical interpretation, the limiting measure is linked to planar maps with 
stars whose type are the monomial of Vi, V 2 and stars of type AB. Those maps are very close from 
Ising configuration on planar graphs except that two stars of type AB can be linked together. For 
integers {kj)i(z{A,B},i<i<D, define 

I{{kj},r,P) = U{ planar maps with kj stars of color i and degree 2j, 

one star of type P (if P 7 ^ 0) and r stars of type AB 
such that there’s no link between any of the r AP-stars. } 

and its rooted counterpart : 

J'{{kj},r,P) = tj{ rooted planar maps with kj stars of color i and degree 2j, 

one star of type P wich is the root and r stars of type AB 
such that there’s no link between any of the r AP-stars. } 
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There’s a relation between these qnantities similar to © : 


We can now relate these numbers to our limit measure : 

Proposition 5.5 Let be as in\5.4\ then on its radius of convergence, 


ht,c{P) = 


1 


1 — 


deg P 
2 


E n 


1 


-f- 






PI \ (1 — j r\ ^ 

kiefflD i,j 3 ’ J 




r^N 


and 


( 19 ) 


F{t,c) - F{t,0) = 


1 -c 2 




-F 


k) 

1 P 


fcieN,ie{i,2},je{i,i3},r>i ^ 


— fci! \ (1 — c^y I F 


,I{{ki},r,0) 


Proof. 

First we define a projection vr from rooted maps to rooted Ising graph such that if M is a map 
7r(M) is obtained by deleting pairs of AB stars which are glued. We now apply Corollary 15.4L and 
translate its result in term of rooted diagrams using m ■■ 


r^N 

All the maps M appearing in that sum are such that tt(M) is an Ising graph rooted at a star of 
type P. For a fixed Ising graph G we must find the contribution in that sum of 7r^“^^(G). But we 
can construct every graph in that set by adding pairs of stars AB on the edges of G. The numbers 
of edges of G is cg = ^ 2jkj so that to get the whole contribution of 7r^“^^(G) we have 

to multiply the contribution of G by 


E 




ai,--- ,aeg,eN 





In that sum, Oj stands for the number of pairs of AB stars added on the i^^ edge. Summing on 
every graphs, we obtain : 





E n 

i,j 

r^N 


-2jt] 

(1 - c2)2i 


\^J{{ki},r,P) 


and the result follows by using m- 

The second point can be proven by proceeding in the same way. 


□ 
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In the rest of this section, we compare a few different approaches to solve the enumeration 
problem of the Ising model. In short, let us emphasize that, for the time being, combinatorial and 
orthogonal polynomials approaches give the more complete and explicit results. However, these 
techniques are still limited to very few models. The Schwinger-Dyson’s equation or the large devia¬ 
tion approaches can be developed for a much wider range of models (such as g-Potts, induced QCD 
etc). However, it seems to us that these arguments still need some mathematical efforts to provide 
as transparent and powerful results (namely for the first a mathematical study of the so-called 
master-loop equations, and for the second a clear understanding of the relations between complex 
Burgers equations and the master-loop equations). A striking difference between the combinatorial 
and the matrix model approaches seems to reside in the fact that matrix models provide for free 
information on the structure of the generating function of the number of planar maps, for instance 
as the Stieljes transform of a probability measure with connected support. 


5.2.1 Orthogonal polynomial approach 

Here we take Vi = V 2 = (g/4)x^. By using orthogonal polynomials techniques, it was proved 
by Mehta |23] that the corresponding free energy Fg^c satisfies 

^g,c - Fo,c = ^ (1 - x)[log f{x) - log 
with f{x) = fg^c solution to the algebraic equation 

/(x){(l - 6-/(x))“^ - c^} -h I2g‘^f{x) - icx = 0 
c 2 

and the root to be taken equals 2“^cx(l — c^)“^ when <7 = 0. 

Starting from there, a simpler expression as been derived in |nj (equation (16), (17) with h = 

z/g) ■ 


Fz.r. = 


-lnh{z) + — 


z-1 


2{3z - 1)3 


+ F 


z + 1 
3z - 1 


j(3z^-3z^F1) 


-h(z)( 


1 




with 


h(z) 


(l-3zr 

1 - c2(l - 3 z)2(1 - 3z2) 


( 20 ) 


Hence, by the preceding, Mehta’s result gives a formula for the generating function of 77 in the 
quadrangulation case. However, it does not a priori gives the limiting spectral measures of the 
matrices. Moreover, this strategy could be only developed completely and rigorously for the Ising 
model and the matrix coupled in chain model m- 


5.3 Direct combinatorial approach 

We can also relate this result to the work of Bousquet-Melou and Schaeffer [7] . Their approach 
is purely combinatorial; they use bijection with well labeled trees (whose generating functions are 
well understood) to obtain algebraic equations for the generating functions of the Ising model. Let 
I(X, V, u) be the generating function of the Ising model on quasi-tetravalent graphs, (i.e. tetravalent 
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graphs except for the root which is bivalent and black) where X (resp. Y) counts the black (resp. 
white) tetravalent stars and u the bicolored edges : 

i(x,r.u)= quasi -tetravalent maps with m tetravalent black stars, 

n tetravalent white stars and r bi-colored edges}. 


If P{x, y, u) is the solution to the algebraic equation : 


P = l + 3xyP^ + 


P{1 + 3xP)(l 3yP) 
u‘^{l — 9xyP^)2 


( 21 ) 


Then, by [7j, Proposition 1 p.4, I can be written in function of P{x,y,u) with x = X{u — and 
y = y(tt - i)2 as 


I{X,Y,u) = + 


P{1 — 3xP — 2xP^ — 6xyP^) yu ^P^(l-|-3xP)^ 


1 — 9xyP‘^ 


(1 - 9xyP2)3 


On the other hand, according to ProDOsition l5.51 if P = tA^+uB^ — cAB and iJ,t,u,c is the associated 
limit measure then on its domain of convergence, 


If we make the following change of variable in CD : 

then we hnd (Unj. Hence, a combinatorial approach can be developed to solve the problem of 
the enumeration of planar maps of the Ising model, a strategy which requires some combinatorial 
insight. The next approach we present, developed in particular by Staudacher, Kazakov and Eynard, 
is a direct analysis of the SD[V] equations. It is a purely analytical and rather robust strategy. 


5.4 Direct study of the SD[P 7 sj„g] equations 

Here, the analysis is based on Theorem 13.41 which asserts that if Pi,P 2 are convex, for small 
parameters, ^ converges almost surely towards the solution of SD[V}^] which is a generating 
function for the enumeration of maps. Hereafter we take c = 1 up to a rescaling x = y/cx, y = y^x, 
Pi(x) = Pi(x), ijj{P{A,X 2 )) = fJ,Jl{P{^/c~^A,^~^X 2 )). Following Eynard [13, we shall analyze 
the solutions of the Schwinger Dyson’s equation. Observe that the following considerations hold for 
any range of parameters, not only small parameters. For large parameters, we do not know that 
the Schwinger Dyson’s equation has a unique solution but we still know that any limit point of the 
empirical measure of the random matrices still satisfies it. In the next section, we shall see that 
for the Ising model and any range of parameters, there is a unique such limit point, and it will 
therefore enjoy the properties described below. We here summarize the main result, as found in 
Eynard m Let Hi be a solution of SD[V/sj„g] 

f,jiiWi{A) - B)P) = f^jiDAP), 
yj{{W'{B) - A)P) = yj0y,j{DBP), 
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with (resp. Db) the non-commutative derivative with respect to A (resp. B) (resp. ^b) and 
Wi{z) = 2:^/2 + Vi{z). Now, let (resp. ^ ib ) be the spectral measure of the matrix A (resp. B) 
then we shall obtain an algebraic equation for H^a{x) (resp. H^b{x)) the Stieljes transform of the 
limiting measure ha (resp. hb) given, for x G C\M by ; 

1 f I 

Hha{x) = Hti - 7 ) = / - duAiv) 

X — A J X — y 

Property 5.6 Let for x^y G C\M, Y{x) = Hij,a{x) — W[{x) and X{y) = Hhb{x) — W 2 {x). Then, 
there exists a polynomial funetion 

d-l 

E{x,y) = ^ aij{t)x^y^ 

so that for all x,y £ C\M 

E{X{y),y)=0 E{x,Y{x)) = 0. 

In particular, ha and hb are absolutely continuous with respect to Lebesgue measure, with Hilbert 
transform Hha and Hhb so that Y{x) = Hha{x) — W[{x) satisfies the same algebraic equation 
with X G M. 

Proof. 

Note that since we know that Ht is compactly supported, we can take in SD[V/sm 3 ] Stieljes 
functions instead of polynomials P since the latest are dense by Weirstrass theorem. 

We take P = P{A) = {x — A)~^ in the second equation in to obtain ; 

=-^ + ^Hta{x) 


Then we use this in the first equation written with 


P{A,B) 


1 [Wl^jy) - Wl^jB)) 
(x -A) {y- B) 


to get after some calculation 


U{x,y){y - Y{x)) = (F(x) - W((x))(x - W'( 2 /)) + 1 - Q{x,y) (22) 


where 


and 


U{x,y) 


( 1 W!2{y)-W!2{B) \ 

^‘\(x-A) (tj-B) )• 


p,,. ,, - M/;(A) W!,(y) - ir'(B)\ 

(i-yi) (y-B) )■ 

To obtain our algebraic equation, we simply define 

E{x,y) = (y(x) - lT((x))(x - W'(y)) + 1 - Q{x,y) 


and we obtain the famous “Master-loop equation” 


E{x, Y (x)) = 0 
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by taking y = Y{x) in (1221) . In a symmetric way, we can show that if X{y) = — W^iy) then 

we also have E(X{y),y) = 0. Note that E is a polynomial function. Hence, this shows that Y(x), 
X{y) and so the generating functions Hy^ix) and HyB{y) are solution to an algebraic equation. 
However, this equation still contains a certain numbers of unknown; {yj{APB'^),p < deg(Vi) —2, q < 
deg(V 2 ) ~ 2}. It is argued in physics that when t is small, the supports of yA and ys should be 
connected and therefore {x,Y{x)) and {X{y),y) should then be genus zero curves. Then, these 
unknowns should be determined by the asymptotic behaviour of X{y) and Y{x) at infinity 

X{y) ~ W 2 iy) - -(1 + o(l)), Y{x) ~ W[{x) - -(1 + o(l)). 
y X 

Note in passing that, as solutions of an algebraic equation, HyA and Hysix) extends con¬ 
tinuously (but in general not differentially) to the real line (eventually as an extended complex 
number). As a consequence, yA and ys have densities with respect to the Lebesgue measure, as 
the limits of the imaginary part of the Stieljes transform on the real line. 

□ 


5.5 Large deviations approach 

A large deviation approach was developed in m, see also Matytsin [22] . Again, we take c = 1 
up to rescaling and denote Wi{x) = x^/2 -|- Vi{x) for z = 1, 2. The main advantage of this strategy 
is to be valid in the whole range of the parameters. Otherwise, it should provide the same type of 
information than in the previous paragraph. Namely, 

Property 5.7 For any polynomials Vi,V 2 going to infinity faster than x^, P-ab converges almost 
surely towards yn = yj which is uniquely defined by the Schwinger-Dyson’s equations 

yj^y-fiDAP) = y-t{{W[{A) - B)P), yj® y-fiOBP) = yAiW^ - A)P) (23) 

and by the fact that yjlA and yj\B (which are the limits of y^ and y^ respectively) are the unique 
minimizers of 

= yA{Wi) + yB{W 2 ) - 2~^ j j \og\x - y\dy'^{x)dy^{y) 

-2“^y j \og\x-y\dy^{x)dy^{y) + ]^m^{j j dxdt + y ^ j pfixfdxdt} 

where the inf is taken over m, p so that yfidx) = pt{x)dx G PO^), yo{x G .) = yA{x G .), 
yi{x G .) = yB{x G .), and 

dtpt{x) + d^mfix) = 0. 

The infimum in {p^,m ) is taken along the solution to a complex Burgers equation; let H = {x G 
M, t G (0,1) : pt{x) > 0} and define on H ut{x) = pt{x)~^mt{x) and ft{x) = ufix) + iirptix). Then 
on Q, 

dtft{x) + ft{x)da:ftix) = 0 . 

Moreover, with yA = y-[\A and yB = yt\B, for yA-almost all x 

W[{x) - uo{x) = HyA{x), yAa.s., W 2 {x) + ui{x) = H yB^x), yB a.s. (24) 

In comparison with the previous statements, we note that the above results hold for all c and Vi, V 2 , 
and not only for small parameters. 
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Proof. 

Most of the proof is contained in m where the convergence of towards the unique 

minimizers of was proved (see Theorem 3.3 in m), as well as the fact that the limit is 

compactly supported and that fij satisfies (1^ but for P G C^(IR) (see section 3.2.1, p. 555 and 
558, in ^n])- It clearly extends to polynomial functions since fij is compactly supported as its 
marginals are. The only point we stress here is that this imply that fij is also uniquely determined. 
Indeed, by proceeding by induction over the degree in i? of a monomial function P, we see that 

t{BP) = -r ® t{DaP) + T{Wi{A)P) 

defines uniquely all the moments t{P{A, B)) from those of t{Q{A)). Note here that this is specific 
to the interaction under consideration ; in general the solutions of SD[V] is not determined by their 
restriction to one variable. 

□ 

Using for instance the fact that if we let gt{x) = tft{x) + x, the Wronskian of (/, g) is null, we 
find that on each connected component of O, there exists an analytic function F so that 

tft{x) + x = F{ft{x)). 

In a small parameter region, it should easily be arguable that is connected, as it is when the 
parameters are null (where the solution at time t can be seen to be a semi-circular variable with 
variance I — t + t'^). According to the previous section, we know that ft extends continuously to 
t = 0 and t = 1 since fiA and hb have densities which yields 

x = F{Mx)) My) + y = F{My)) (25) 

for all X in the support of jjLA and all y in the support of /r^. Noting that fo{x) = W{{x) — HyA{x) = 
—Y (x), fi{x) = —WJ^^ix) +HyiB{x) = —X{x) it is tempting to hope that (|^ yields the same result 
that Property on namely that {Y{x),x) and {y,X{y)) satisfy the same algebraic equation. Our 
knowledge of this field is much too limited to unable us to get this conclusion. 
Acknowledgments : We are extremely grateful to B. Eynard and G. Schaeffer for many comments 
which helped us to compare the different mathematical approaches to the enumeration of planar 
maps. We also thank A. Okounkov for many useful discussions. 
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